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Abstract. 

We discuss transport through double gated single and few layer graphene devices. 
This kind of device configuration has been used to investigate the modulation of the 
energy band structure through the application of an external perpendicular electric 
field, a unique property of few layer graphene systems. Here we discuss technological 
details that are important for the fabrication of top gated structures, based on electron- 
gun evaporation of Si02- We perform a statistical study that demonstrates how - 
contrary to expectations- the breakdown field of electron-gun evaporated thin Si02 
films is comparable to that of thermally grown oxide layers. We find that a high 
breakdown field can be achieved in evaporated Si02 only if the oxide deposition is 
directly followed by the metallization of the top electrodes, without exposure to air of 
the Si02 layer. 
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1. Introduction 

The ability to isolate and embed single- and multi-layer graphene in double gated 
structures is paving the way to reveal unique electronic properties of these systems 
[HElEllllISlElITlIHlini HUl HH H2l H3]. The ability to change the voltages applied 
to a nano-fabricated top gate and to the back-gate offers the possibility to gain local 
control of the charge density and of imposing locally a perpendicular electric field. This 
device configuration was used recently to show how the band structure of graphene-based 
materials can be tuned continuously [8], [T2] . In particular, bilayer graphene exhibits an 
electric field induced insulating state due to the opening of a gap between valence and 
conduction band [8], and in trilayers, which are semi- metals, the band-overlap can be 
increased substantially [12] . Cleverly designed top gates on a graphene single layer have 
also been used successfully for engineering p-n junctions [H EJ [31 El EJ [7] , necessary for 
the investigation of Klein tunneling [HOE], and to attempt the fabrication of controlled 
quantum dots [151 [IB]- Key to the fabrication of top gated structures is the ability to 
deposit good quality thin gate oxides, with high breakdown field and low leakage current. 
Here, after reviewing the relevance of double gated devices on the electric field 
modulation of the band structure of double and triple layer graphene, we discuss in some 
details technological aspects related to the properties of the Si0 2 layers used as gate 
insulators. In particular, we discuss how we can routinely achieve high breakdown fields 
in electron-gun evaporated thin SiO"2 films (15 nm), comparable to the breakdown fields 
of thermally grown SiO"2, which is surprising given that SiC>2 deposited by evaporation 
was long believed to be a poor quality insulator. To unveil the reasons behind the 
good insulating quality of our evaporated SiO"2 films, we conducted a statistical study 
of leakage current and breakdown voltage in capacitors, where two metallic electrodes 
are separated by a SiO"2 layer fabricated in different ways. We demonstrate that if Si02 
and top gate metal electrodes are deposited subsequently without exposing the SiC>2 to 
air, the electrical performance of electron-gun evaporated Si02 is comparable to that of 
thermally grown SiC>2. In contrast, exposure to air of the Si02 layer before deposition 
of the counter-electrode leads to much worse insulating characteristics. Our findings 
indicate that extrinsic degradation -probably due to the absorption of humidity- has 
limited in the past the insulating quality of electron-gun evaporated Si02- 

2. Device and fabrication 

Single and few layer graphene flakes used for the device fabrication were obtained 
by micro-mechanical cleavage of natural graphite crystals, and by their subsequent 
transfer onto a highly doped Silicon substrate (acting as a gate) covered by a 285 nm 
thick thermally grown Si02 layer. The thickness of the graphene layers can be simply 
identified by analyzing the shift in intensity in the RGB green channel relative to the 
substrate (i.e. Relative Green Shift) [H [121 HH EE3 EH] • A plot °f the relative green shift, 
as extracted from optical microscope images of various samples taken with a digital 
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camera, exhibits plateaus corresponding to the discrete thickness values -see Fig. [Th. 
Subsequent transport measurements (quantum Hall effect, resistance dependence on a 
perpendicular electric field, etc.) confirm the validity of this optical method (Fig. [lb). 
The fabrication of nanostructures is accomplished by conventional electron-beam 
lithography. Metallic contacts and top gates were deposited by electron-gun evaporation, 
respectively of Ti/Au (10/25nm thick with a back ground pressure of 9 * 10 _7 torr) for 
contacts and SiC^/Ti/Au (15/10/25nm thick with a back ground pressure of 3 * 10 -7 
torr) for top gates, followed by lift-off. We took special care to fabricate all the ohmic 
contacts within 60 nm from the edges of the top gated areas, so that two probe resistance 
measurements are dominated by the resistance of the double gated region [19]. All 
transport measurements in double gated devices (see Fig. [lb) were made using a lock-in 
technique (excitation frequency: 19.3Hz), in the linear transport regime, at temperatures 
ranging from 300 mK up to 150K. 

To understand why, contrary to expectations, we manage to achieve high breakdown 
field in thin, electron-beam evaporated SiC>2 films, we conducted a macroscopic study of 
the breakdown characteristics on two types of capacitor test structures. The first -which 
we refer to as type A- is characterized by subsequent evaporation of SiC^/Ti/Au without 
breaking the vacuum in between the deposition of the different materials. For the second 
-type B Si0 2 - we exposed the device to ambient for one hour after the Si0 2 deposition, 
before evaporating the Ti/Au counter electrode. The breakdown test measurements 
were made with a Keithley-2400 source-meter on more than 130 different capacitors 
(with three different surface areas: 125 x 115/xm 2 , 175 x 150/im 2 and 215 x 195/im 2 ). 

3. Transport experiments in double gated few layer graphene devices 

The measurement configuration used for double gated devices is shown in Fig. [lb. A 
finite voltage applied to either one of the gates (back or top gate) changes the position 
of the Fermi level in the corresponding gated region of the graphene layer, by an amount 
corresponding to the induced charge density. In addition, by biasing the two gate elec- 
trodes with opposite polarity, a large external electric field applied perpendicular to the 
layer is generated, which is equal to E ex = (V bg — V tg )/d to t (d to t = 15 + 285nm is the 
total SiC>2 thickness). In this device configuration, we can monitor the evolution of the 
in plane transport properties for each few layer graphene device as a function of E ex . 
Fig. [Tbl, e and f show the typical behavior of the square resistance measured respec- 
tively in double gated single (d), double (e) and triple (f) layer graphene, when sweeping 
the top gate, while keeping the back gate at a fixed potential. It is apparent that the 
overall electric field dependence of R sq is markedly distinct for graphene layers of differ- 
ent thickness. In all cases, the resistance exhibits a maximum {R™ ax ) whose value and 
position in gate voltage depend on the voltage applied to the gate on which a fixed po- 
tential is applied during the measurement. At E ex = V/m we find that for single and 
bilayer graphene R™q X ~ 6KQ, close to a conductance per square of 4e 2 //i, as expected, 
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Figure 1. Panel (a) shows a plot of Relative Green Shift (as denned in the main 
text) for 90 few layer graphcnc areas, showing clear plateaus for graphcne layer of 
different thickness. Panel (b) shows the conductance of a device (indicated by the 
blue dot in the plot shown in panel (a)) measured in the presence of high magnetic 
field (B=8T); plateaus are observed at conductance values characteristic for Dirac 
electrons, indicating that flakes with Relative Green Shift = 0.06 are indeed single 
layers. Panel (c) shows a schematic representation of double-gated graphcnc devices, 
together with the measurement configuration used to study transport as a function of 
voltages applied to the top and back gates. Panels (d), (e) and (f) show plots of the 
square resistance in single (d), double (c) and triple (f) layers respectively, measured 
while sweeping the top gate, for different fixed voltages applied on the back gate (for 
single layer T = 300mK, W = 2ftm, L = lAftm, ft = 3500cm 2 /Vs, double layer 
T = 300mK, W = 1.7/xm, L = lfjtm, ft = 900cm 2 /^s and for triple layer T = 1.5K, 
W = lftm, L = lAftm, ft = 800cm 2 /Vs). 



indicating that the fabrication of top gate structures does not damage significantly the 
material (for trilayers the square resistance is somewhat lower, owing to the presence 
of an overlap between valence and conduction band). Increasing the external electric 
field induces a well defined -and different- response for the square resistance of layers of 
different thickness. Respectively, in a single layer -R™ aa: is not affected by a finite E ex ; in 
bilayers at low temperature i?™ ax increases from 6KQ to very large values (> 100 Kf2); 
in trilayers i?™ ax decreases with increasing E ex . These experimental findings, confirmed 
in a number of different samples (3 single layers, over 10 double and 10 triple layers) 
provide a clear indication that each few layer graphene is a unique material system, with 
distinct electronic properties. 

Transport measurements over a wide range of temperatures (from 300mi^ up to 150K) 
underline the unique electronic properties of these few layer graphene devices. In bilayer 
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graphene the larger E ex , the more pronounced is the temperature dependence of -R™ ax , 
see Fig. [2^ and b. At E ex = V/m, R™ ax is only weakly temperature dependent (as it 
is typical of zero-gap semiconductors), and at E ex ^ V/m the observed behavior is the 
one typical of an insulating state. On the contrary, trilayer graphene devices display a 
decrease of R™ q ax when lowering the temperature, stemming for the semimetallic nature 
of the constituent material. 

At a more quantitative level we find that R^ x in bilayer graphene is well described 
by oc exp(T /T) 1 / 3 , with T « 20 K at the maximum applied external electric field (see 
inset in Fig. [2b). This temperature dependence is indicative of variable-range hopping 
in a two dimensional material where an energy gap has opened, and where disorder 
causes the presence of sub-gap states (T is related to this subgap density of states [8]) 
-making it difficult to estimate A from transport experiments. 

Indeed, in a disorder free bilayer graphene, at E ex ^ a gap (A) opens in the band struc- 
ture and the density of states in the gapped region is zero. In this ideal case, when the 
Fermi level is in the middle of the gapped region the value of R™ ax at a finite temperature 
is entirely determined by thermally activated charge carriers {R™ ax oc exp(A/2/c j e7 1 )). 
Therefore, A can be accurately determined from a plot of Zn(i?™ ax ) versus 1/T. On the 
other hand, in real devices the presence of disorder creates a finite density of states in 
the the band gap of bilayer graphene. Now charge carriers can conduct via variable- 
range hopping at -R™ ax . In this case R™ ax is a function of the density of states at the 
Fermi level and not any more simply a function of A. 

The fact that a gap opens in the band structure of bilayer graphene at finite E ex is 
evident from the temperature dependence of R sq as a function of charge density. In 
particular, when the Fermi level lays deep into the conduction and/or valence band, 
a temperature independent R sq is expected. However, when the Fermi level is shifted 
across the energy gap region R sq should display an insulating temperature dependent 
behavior. Experimentally this is achieved by measuring R sq at a fixed value of either 
of the gates (e.g. Vb g ) and for different voltage applied on the other gate (e.g. V tg ) (see 
Fig[2|i and b). The cross over from band transport to variable hopping range in the gap 
occurs at the edge of the valence and conduction band (see Fig. |2b). 
Similar previous transport experiments in double gated bilayer [8] reported an energy 
scale of l-10mV associated with the insulating state induced by E ex ^ 0. This energy 
scale seen in transport is much smaller than the energy gap recently probed in opti- 
cal spectroscopy experiments (A ~ 200mV at E ex = 2V/nm [22]). Possibly the finite 
sub-gap density of states induced by the disorder is at the origin of the small energy 
scale measured in transport experiments, however the specific mechanism responsible 
for these experimental observations remains an open question. 

The temperature dependence of the resistance in trilayer graphene is opposite to the one 
observed in bilayers, and it reveals that this material system is a semimetal with a finite 
overlap (5e) between conduction and valence band. This band overlap can be estimated 
within a two band model [121 EDI El] , where the number of thermally excited carriers 
increases with temperature according to n(T) = (IQnm* /h 2 c)kBTln(l + e fe / 2i ^s T ) [m* 
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Figure 2. Panels (a) and (c) show plots of the square resistance at different 
temperatures for the bilayer and trilayer devices whose data are shown in panel (e) 
and (f) of Fig[TJ Panel (b) shows a plot of R sq vs. T and ln(R sq ) vs. T 1 / 3 extracted 
from the measurements shown in panel (a), at Vb g = — 50V and different Vt g values 
as indicated in the legend. The external electric fields applied on the bilayer for each 
different V tg are respectively -0.177 V/nm (■), -0.175 V/nm (;. ), -0.173V/nm (•), 
-0.167 V/nm (a), -0.16 V/nm (t), -0.153 V/nm (■*). Panel (d) shows the normalized 
charge density as a function of temperature, extracted from the trilayer measurements 
of panel (c). The continuous line is a fit based on a two band model with finite overlap 
(see main text). 



effective mass and c equal to twice the layer spacing). Measurements at finite E ex show 
that 5e decreases when increasing external electric field (5e goes from 32 meV to 52 
meV in the measurements of Fig. [2H). 

These experiments demonstrate that a perpendicular electric field applied on few layer 
graphene is a valuable tool to change the band structure of these materials. Double 
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gated structures lead to the discovery of the only known electric field tunable insulator, 
i.e. bilayer graphene, and of the only known electric field tunable semimetal, i.e. trilayer 
graphene. 

4. Evaporated Silicon oxide as top gate dielectric 

The opening of a sizeable band gap in bilayer graphene, and large changes in the band 
overlap of graphene trilayers require the application of large external electric fields to 
these material systems. It is the breakdown field of the gate-dielectric that imposes 
a limit to the maximum value of E ex experimentally accessible in double gated struc- 
tures. To optimize this aspect of the devices, we conducted a systematic study of the 
breakdown electric field of SiC>2 gate oxide for devices with different areas, fabricated 
under different conditions. Here we discuss the details of this investigation. From our 
statistical analysis we conclude that what is crucial is not the SiC>2 deposition method, 
but the details in the metallization of top gate electrodes afterwards, which affect the 
final quality of the oxide gate dielectric. 

We compare capacitor devices fabricated following two different procedures for one of 
the oxide dielectric/metal electrode interfaces. In particular, the devices were fabricated 
on a Si/Si0 2 substrate (identical to the one used for the graphene devices previously 
described) on which we deposit a Ti/Au (5/20 nm) film -common electrode for the ca- 
pacitors. In devices of type A, the Si0 2 deposition and top electrode (Ti/Au 5/20 nm) 
metallization processes were carried without breaking the vacuum. On the contrary in 
devices of type B, the SiC>2 gate dielectric was exposed to air for one hour prior to the 
deposition of the top electrode metals. The SiC>2 deposition was carried out typically 
at 3 * 10~ 7 torr back ground pressures. We did not observe a dependence of the insul- 
tating properties of the Si0 2 dielectric, breakdown field and leakage current over the 
background pressure range from 1 * 10~ 7 to 5 * 10~ 7 torr. 

Fig. [3^ and b show various I — V traces, measured in ambient condition, for type A and 
type B devices. A first clear difference between the two types of devices is the magnitude 
of the leakage current, visible by plotting the I — V curves both in linear and logarithmic 
scale, see Fig. [31 For a surface area of 215 x 195yum 2 , we find leakage — 2-5 * 10~ 10 A/[im 2 
for the best type B devices which is one order of magnitude larger than that measured 
in the worst type A devices (the difference for typical devices are much larger than one 
order of magnitude). This extremely different level of leakage current already indicates 
that the exposure of SiC>2 to air previous to the deposition of top metals has a large 
negative influence on the insulating performance of the oxide. 

The I — V characteristics further show that for a fixed surface area (215 x 195/im 2 ), the 
breakdown voltage Vbd f° r type A is typically in the range 8V < Vbd < whereas 
type B devices breakdown anywhere in the range 0V" < Vbd < 6 V. The differences 
in the failure of device types are best summarized in the histogram plots of Vbd -see 
Fig. Hh- For type B we find a large spread in the distribution of Vbd, in contrast to 
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Figure 3. Panels (a) and (b) show I — V characteristics of respectively type B devices 
with S = 215 x 195/um 2 (a), and type A devices with S = 215 x 195//m 2 (blue curves), 
S = 175 x 150/xm 2 (red curves) and S = 125 x 115/zm 2 (black curves) (b). Panels (c) 
and (d) show the I — V curves in logarithmic scale for the two types of devices and 
areas S = 215 x 195/^to 2 (black curves in (c) (d)) and S = 175 x 150fim 2 (red curves 
in (d)). 

the narrow distribution characteristic of type A devices. Furthermore, the comparison 
of Vbd for type A devices with different surface areas shows that Vbd increases slightly 
with decreasing the device area, possibly indicating that the properties of SiC>2 in type A 
close to breakdown are determined by small defects present in the film with rather small 
probability. However, we cannot rule out that the differences between the different sam- 
ple populations originate from small differences (2-3 nm) in the thickness of the SiC>2 
layers. Note, in fact, that the leakage currents of these devices at low bias have only 
small sample-to-sample fluctuations, suggesting that the SiC>2 layers in type A devices 
are very uniform. 

To try to quantify better our observations and analyze the role played by the spe- 
cific fabrication technique and surface areas on the device performance (e.g. breakdown 
field), we adopt a failure analysis methodology [24]. In what follows we provide a sta- 
tistical description of the breakdown probability introducing the cumulative probability 
(P) as the probability of a device to breakdown at a given voltage. From failure method- 
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Figure 4. Panel (a) shows histogram plots of the breakdown field for different 
populations of test devices as specified in the legends. The panels in (b) and (c) 
show the cumulative probability for type A devices for different surface areas (dots are 
experimental data, and the continuous line is a fit to the Weibull distribution). The 
plot in d) is a fit to the cumulative probability for type B devices (S = 215 x 195^m 2 ) 
to the Weibull distribution. 

ology, we notice that possibly the most flexible distribution for the failure of a population 
of samples, is the Weibull distribution [231 125] P = 1 — exp[— {Vbd /VqY S / Sq] (S capac- 
itor surface area, S reference surface area, (3 and V are the Weibull parameters). The 
parameter j3, also known as Weibull shape parameter, determines the shape of the prob- 
ability density function -i.e. higher (3 indicates distributions with low dispersion of Vbd- 
Vq is the Weibull scale parameter, whose only effect is to scale the Vbd distribution (the 
larger Vo, the more "stretched" the distribution). Depending on the value of the Weibull 
parameters, this distribution mimics the behavior of other statistical distributions such 
as the normal and the exponential. Given a sample population, the Weibull parameters 
provide a quantitative measure of the failure probability. Both (3 and Vo are strongly 
affected by the failure mechanism which can eventually be identified when comparing 
the Weibull parameters for different sample populations. For instance, a value of (3 that 
does not depend on the capacitor surface area -i.e., the variance does not change with 
the surface area-, means that the microscopic mechanism of breakdown is common to 
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all the samples, independent of the specific area [25J. 

The good agreement between a fit to the Weibull distribution of the cumulative prob- 
ability for each different surface area and device types shows that the data of each 
different device population is well described by the Weibull distribution. To evaluate 
whether a single Weibull scaling law can explain breakdown results for all the different 
surfaces, we notice that -for type A- we can fit all the cumulative probability distri- 
butions with the same value (3 = 52, see Fig. 0b. This is made apparent by plotting 
ln(— ln(l — P)) = ln(S/So) — flln(Vo) + (31ti(Vbd) as a function of Iti(Vbd) for data 
sets corresponding to different surface areas. Fig. [4b and c show clearly that all the 
data are lined up on parallel linear slopes (i.e., (3 is the same in all these cases). This 
finding implies that the breakdown mechanism is the same for all studied surface areas. 
Note that Vq varies slightly with surface area, and that, as mentioned above, we cannot 
exclude that the origin of these variations is a small difference in the thickness of the 
Si0 2 layers for the different device populations (a difference of 2-3 nm would suffice to 
explain this observation). 

We notice that /5 = 52 -estimated from the fit in Fig. Hb and c- is comparable to values 
found for thermally grown thin SiC>2 films, and it is compatible with failure of the de- 
vices due to surface roughness [23] probably being transferred to the dielectric film from 
the Ti/Au substrate. This quantitative analysis make it possible to state that electron- 
beam evaporated SiC>2, directly coated by a metallic layer without exposure to air has 
an essentially identical quality to that of thermally grown oxide. A similar analysis of P 
for type B devices, gives a (3 = 3.3 -i.e. a much higher dispersion of breakdown field (see 
Fig. HU). This small value for (3 in type B devices quantifies the much larger statistical 
spread of the oxide properties in these devices. Since the only difference between type A 
and B devices is the fabrication step of the SiC^/Ti interface, we conclude that exposure 
of the SiC>2 to air is indeed the cause for the poor insulating qualities. Indeed it is well 
known that SiC>2 is an hygroscopic material, that easily absorbs humidity in air. The 
humidity absorbed can affect the composition of the entire layer providing paths for 
the leakage current, and creating weak spots at which breakdown occurs already at low 
voltage. 

5. Conclusions 

In conclusion we have briefly reviewed transport in double gated bilayer and trilayer 
graphene devices. Motivated by the need for large electric fields, we have conducted a 
statistical study of the breakdown field for over 100 top gated structures fabricated in 
different conditions and with different surface areas. Adopting a failure analysis based 
on the Weibull distribution, we show that the most reliable top gates are obtained when 
depositing in SiC^/Ti/Au without breaking the vacuum. Electron-beam SiC>2 layers 
evaporated in these ways have insulating characteristics as good as those of thermally 
grown SiC>2 layers. 
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